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ABSTRACT 

Mass loss of red supergiants (RSG) is important for the evolution of massive stars, but is not fully explained. Several 
empirical prescriptions have been proposed, trying to express the mass-loss rate (M) as a function of fundamental stellar 
parameters (mass, luminosity, effective temperature). Our goal is to test whether the de Jager et al. (1988) prescription, 
used in some stellar evolution models, is still valid in view of more recent mass-loss determinations. By considering 40 
Galactic RSGs presenting an infrared excess and an IRAS 60-fim flux larger than 2 Jy, and assuming a gas-to-dust 
mass ratio of 200, it is found that that the de Jager rate agrees within a factor 4 with most M estimates based on the 
60-/im signal. It is also in agreement with 6 of the only 8 Galactic RSGs for which M can be measured more directly 
through observations of the circumstellar gas. The two objects that do not follow the de Jager prescription (by an 
order of magnitude) are fj, Cep and NML Cyg. We have also considered the RSGs of the Magellanic Clouds. Thanks 
to the works of Groenewegen et al. (2009) and Bonanos et al. (2010), we find that the RSGs of the Small Magellanic 
Cloud have mass-loss rates consistent with the de Jager rate scaled by (Z/Zq) 01 , where Z is the metallicity and a is 
0.7. The situation is less clear for the RSGs of the Large Magellanic Cloud. In particular, for L > 1.6 x 10 5 Lq, one 
finds numerous RSGs (except WOH-G64) having M significantly smaller than the de Jager rate and indicating that 
M would no longer increase with L. Before this odd situation is confirmed through further analysis of LMC RSGs, we 
suggest to keep the de Jager prescription unchanged at solar metallicity in the stellar evolutionary models and to apply 
a (Z/Zq) ' 7 dependence. 
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1. Introduction 

The red supergiant stars (hereafter RSG) are massive stars 
evolving beyond the main sequence and passing a fraction 
of their lives in the cool upper region of the Hcrtzsprung- 
Russell (H-R) luminosity-temperature diagram. In the so- 
lar neighborhood, the initial mass of RSGs is about 8 
to 25 M , their effective temperature is for most objects 
3000 to 4000K, and their luminosity is between 2 x 10 4 
and ~ 6 x 10 5 L Q . Several important works have re- 
cently revisited the basic properties of these stars in the 
Galaxy and in other systems like the Magellanic Clouds and 
M31 (Massey 1998; Levesque et al. 2005, hereafter LM05; 
Massey et al. 2008; Levesque 2009 and references therein). 
In general, a rather good agreement between recent obser- 
vations and the predictions of evolutionary models has been 
established. However, one of the important processes occur- 
ing at the surface of these stars, mass loss, is still poorly 
known. The radii of RSGs are huge, typically 500-1500 R©. 
Consequently, the gravity of these stars is very small, about 
50 000 times smaller than that of the Sun. The small grav- 
ity favours ejection of matter to the interstellar medium, 
and indeed some RSGs lose mass at a high rate, 10~ 6 to 
10~ 4 Moyr- 1 . 

Mass loss has several consequences. First, a significant 
fraction of the stellar mass is lost, which is important for 
subsequent evolution of the star. Secondly, the ejected cir- 
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cumstellar matter enriches the ISM in dust, e.g. diverse 
types of silicates (see for example Verhoelst et al. 2009) 
and even carbonaceous matter (Sylvester et al. 1994, 1998). 
One can note also that some types of core-collapse super- 
novae have RSG progenitors (Smartt et al. 2009; Josselin & 
Langon2009), so that the density and radial distribution of 
the RSG winds play an important role in the expansion and 
evolution of the supernova remnant (Chevalier et al. 2006; 
Smith et al. 2009; van Loon 2010). 

Despite the importance of RSG mass loss, a quantita- 
tive physical model of the winds has yet to be established. 
Like for other mass-losing cool giants, atmospheric motions 
such as pulsation, convection and shocks occur (e.g. Kiss 
et al. 2006; Josselin & Plez 2007), and this may produce 
out-moving gas sufficiently cool for dust to form Then, ra- 
diation pressure on the grains may be sufficient to drive 
the outflow (see for example Habing et al. 1994). Although 
this scenario is plausible, some important details are still 
unclear. For example, the condensation of the Fe-rich sil- 
icate grains, crucial to drive the wind, is not well under- 
stood (Woitke 2007). When dust is incompletely formed 
or absent, other mass-loss mechanisms have to be found, 
but none is completely satisfactory (e.g. Bennett 2010 and 
references therein). 

Because no physical model exists, diverse empirical 
parametrizations are used in stellar evolution codes for esti- 
mating M as a function of basic stellar characteristics, i.e. 
mass, luminosity, and effective temperature. An important 
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example is the empirical law built by de Jager et al. (1988; 
hereafter dJ88). This law was based on observational esti- 
mates of M for stars located over the whole H-R diagram, 
including 15 RSGs. It is used in the Cambridge stellar evo- 
lution code with a (Z/Zq) 5 metallicity scaling (Stancliffc 
& Eldridge 2009). The dJ88 law was also used in past evolu- 
tion models of the Geneva group (described in, e.g., Maeder 
& Meynet 1989; Eggenberger et al. 2008). However, in the 
most recent Geneva code, the dJ88 prescription is used only 
for yellow supergiants with T e g between 5000 K and 8000 K. 
For smaller T c h, a power-law in L, M oc L 17 , is adopted 
with no metallicity scaling. 

The goal of this paper is to investigate whether this 
de Jager formulation is still in agreement with mass-loss 
rates obtained more recently. For example, mass-loss rates 
of RSGs were estimated with IRAS observations of the cir- 
cumstellar dust (Jura & Kleinman 1990, hereafter JK90; 
Josselin et al. 2000, hereafter JB00) or from ISO spectra 
(Verhoelst et al. 2009). Several radio millimeter observa- 
tions of CO from RSGs have also been made since 1988 
(JB00; Woodhams 1993). The question is whether these 
more recent measurements support the dJ88 prescription. 

In the following, Sect. 2 recalls the prescription of de 
Jager and compares it with a few others. In Sect. 3, we 
compare the dJ88 rate with measured mass-loss rates of 
RSGs, and we consider three different cases: i) the RSGs 
of the solar neighbourhood and M estimated with their 
infrared excess at 60/im; ii) the eight cases of these RSGs 
for which circumstellar gas has been directly observed and 
used to derive M ; hi) the RSGs of the Magellanic Clouds. 
In Sect. 4, we discuss our results and conclude in Sect. 5. 

2. Comparison of de Jager rates with other 
prescriptions 

In the paper of dJ88, the authors collected a large num- 
ber of mass-loss rates determined for 271 Galactic stars 
located over the whole H-R diagram. The mass-loss rates 
could be matched by a single empirical formula in which 
the parameters are the effective temperature T g and the 
luminosity L. Hereafter, we call this formula "the de Jager 
rate". According to dJ88, the accuracy of this formula is 
about 0.45 in log M, i.e. about a factor of 3. Because the 
de Jager formula comprises 20 Chebycheff polynomia, the 
dependence of M on T c g and L cannot be seen clearly. 
Therefore, we show in Fig. 1 what is the situation for RSGs. 
The de Jager rate is plotted for two different effective tem- 
peratures, 3500 and 4000 K, since the majority of Galactic 
RSGs have their T e g between these values (LM05). More 
precisely, in the list of 74 RSGs of LM05 (with a M-type 
classification), only a Her has a lower temperature, with 
3450 K, but it is not a massive star, M w 2.1 M Q accord- 
ing to Schroder & Cuntz (2007). Figure 1 shows that M 
is not very dependent on T e g in the de Jager formulation, 
and that M depends only on L. 

Among the 271 stars (hot or cool, luminous or not), 
which are at the basis for the dJ88 prescription, there are 
only 15 Galactic RSGs, for which dJ88 compiled T e g, L and 
an observational estimate of M, or an average of estimates. 
These compiled data are listed in Appendix A, together 
with results of this work, discussed later. It can be noted 
that the T e g are significantly cooler than those recently de- 
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Fig. 1. Comparison of several prescriptions for the RSG 
mass-loss rates (unit M Q yr _1 ) plotted against luminosity 
(in L Q ). The de Jager prescription is represented by the 
2 curved lines (solid line T cS = 4000 K, dotted line T eff = 
3500 K). 

termined by LM05, when available. In addition, about half 
of the estimated M are based on the early paper of Gehrz 
and Woolf (1971). So, we think it appropriate to reexamine 
the validity of the dJ88 prescription with a larger sample 
of objects (40 are considered in this work) and more recent 
data. 

There exists a second formulation published later by 
Nicuwenhuijzen and de Jager (1990; hereafter NJ90). The 
stellar mass is taken into account, and the dependence of 
M on stellar mass, luminosity and temperature is more 
explicit. NJ90 find: 

log M = -7.93 + 1.641og L + 0.16 log M - 1.61 log T cff . 

It can be seen in the above formula that M depends 
very little on M, since for M =9 M and 25 M Q , the second 
term is equal to 0.15 and 0.22 respectively. This represents 
a difference of only 17 percent on M , and consequently this 
mass term can be considered to be constant and equal to 
0.18. Then, a more explicit formulation is: 

M = 5.6 x 10- 6 (L/10 5 ) 164 (T cff /3500)- 1 61 . 

The rate of NJ90 is plotted in Fig. 1 for T cff = 3750K. 
For other temperatures, this rate does not change much: 
for effective temperatures of 3500 or 4000K, the NJ90 line 
in Fig 1 would be shifted up or down by only 0.05 in the y- 
axis. This NJ90 is not very different from the de Jager rate: 
the largest difference is about a factor of 2 at L ~ 10 5 . NJ90 
also note that for cool evolved stars the dJ88 rate is more 
accurate than the NJ90 one in fitting the M observational 
measurements (their Table 2). 

We have also plotted in Fig. 1 the rate adopted in the re- 
cent Geneva models (R. Hirschi, private communication). 
This rate is represented by a continuous straight line lo- 
cated very close to, and just below the NJ90 line. This 
Geneva rate is practically identical to the NJ90 one, but 
there is no dependency on T e g . The Geneva rate is the fol- 
lowing: M = 4.7 x lO-^L/lO 5 ) 1 - 7 . 



2 



N. Mauron and E. Josselin: The mass-loss rates of red supergiants 



Figure 1 also shows the Reimers law (Reimers 1975; 
Kudritzki & Reimers 1978): 

M = 5.5 x 10~ 13 LR/M, 

with L, R and M in solar units. Since we restrict our sur- 
vey to RSGs, this relation can be simplified: R is equal to 
^0.5 x (y off y5770)~ 2 and the factor depending on T e g does 
not vary much (2.72 and 2.08 for T cff = 3500 and 4000 K 
respectively). In addition, for the RSGs studied here, there 
is an approximate relation between mass and luminosity. 
For M = 9M Q , L w 2.5 x 10 4 L Q , and for M = 25 M , 
L w 3 x 10 5 L Q (see, for example, Fig. 8 of Massey et al. 2009, 
lower panels). We find M«0.14L - 41 . So the Reimers rela- 
tion for an average temperature T e g = 3750K is: 

M = 2.6 x 10- 6 (L/10 5 ) 1 - 09 . 

This relation is plotted in Fig. 1, and it can be seen 
that the Reimers rate is close to the dJ88 rate (within a 
factor of 2 for L < 2 x 10 5 L Q ), but the latter is signifi- 
cantly stronger for larger luminosities. Again, we have not 
drawn, for clarity, the Reimers lines corresponding to other 
temperatures, but for 3500 and 4000K, they would be 0.060 
(in log) above and 0.056 below the plotted line. So the de- 
pendance on temperature is quite small compared to the 
observational uncertainty on measured mass-loss rates con- 
sidered in the following (typically a factor 2 or 3). 

There exist other mass-loss prescriptions. Feast (1992) 
noted that there is a very good correlation between the 
pulsation period P of 15 RSGs of the Large Magellanic 
Cloud and their luminosities L based on JHK and IRAS 
data and estimated by Reid et al. (1990). The mass-loss 
rates published by Reid et al. for these stars also correlate 
rather well with P. A relation linking M and L is derivable 
and given in Equation (12) of Salasnich et al (1999). This 
relation (plotted in Fig. 1) can be written: 

M = 2.1 x lO^OL/lO 5 ) 1 - 385 . 

A prescription based also on the LMC RSGs and 
the work of Reid et al. (1990) is used by Vanbeveren et 
al. (1998) (see also Vanbeveren et al. 2007). It is labelled 
"Vanb." in Fig. 1. This rate is in agreement with the Feast 
rate for L ~ 10 5 L Q , but its slope is lower. Its expression is: 

M = 2.0 x lO-^L/lO 5 ) - 80 . 

Another mass-loss prescription has been proposed by 
Salasnich et al. (1999), who adopt the Feast relation, but 
take a gas-to-dust mass ratio increasing with luminosity. 
Typically, this ratio would be 200 for AGB stars but could 
increase to 1000 for RSGs with a luminosity of 10 5 L Q . 
Consequently, the Salasnich prescription implies very high 
mass- loss rates. For L ~ 10 5 L Q , this prescription gives 
10~ 4 M Q yr _1 , about a factor of 50 above the de Jager 
rate. The Salasnich law is: 

M = 1.00 x lO^OL/lO 5 ) 2 - 1 . 

Finally, van Loon et al. (2005) derived a mass-loss law 
for AGB stars and RSGs, but it cannot be applied if the 
circumstellar envelope is not very dusty. Nevertheless, the 
rate of van Loon et al. (2005) is interesting to compare with 
the de Jager rate, and this will be done with more details 
in the discussion. In the next section, we compare the de 
Jager rate with modern mass-loss measurements. 



3. Comparison of the de Jager rate with 
observations 

3.1. The 60-fim excess emission of RSGs in the solar 
neighbourhood 

The RSGs of the solar neighbourhood are key targets for 
investigating mass loss, either because they are close to us 
like Betelgeuse or Antares, or because the membership to 
an OB association gives a distance indication and confirms 
the object to be a massive star. These RSGs have been orig- 
inally listed by Humphreys (1978). JK90 investigated mass 
loss of the most luminous RSGs. They consider a list of 21 
objects more luminous than 10 5 L Q and estimate mass loss 
with the IRAS 60^m flux. JB00 consider 65 RSGs, present 
photometric JHKL measurements of most of them, and de- 
tect the CO millimeter emission in several cases. One of 
their conclusions is a large scatter of M for a given lumi- 
nosity, a property that was also noted by JK90. LM05 have 
revised effective temperatures, distances and interstellar ex- 
tinction for 74 RSGs. 

A first list of RSGs was built from the list of 74 
objects of LM05, supplemented with 21 objects from 
JB00. For these 95 sources, we collected optical (UBVI) 
and near- infrared photometry (JHKL), preferentially from 
Lee (1970). Data from Simbad, 2MASS, JB00 and Gezari 
et al. (1993) were also used. At mid-infrared wavelengths, 
IRAS measurements were adopted. Then, for most of the 
stars, photometry from the U band to the 60^im band is 
available. Almost all these stars belong to OB associations, 
and the distance moduli and interstellar absorption given 
by LM05 were adopted. When a star was not examined in 
the work of LM05, we took the distance modulus of the 
association from Humphreys (1978) and the reddening was 
derived from the averaged reddening of the OB-type stars 
of the association. For several objects (e.g., VY CMa, VX 
Sgr), the distance and absorption come from recent, indi- 
vidual determinations. 

From this primary list of 95 objects, we are interested 
in those presenting evidence of mass loss. For that purpose, 
we selected objects presenting an infrared excess, measured 
by the ratio of 12 /im flux to the dereddened K-band flux. 
Following JK90, only stars with /u/Jk larger than 0.1 
were selected. In addition, we imposed feo > 2 Jy, in or- 
der to work with reliable 60/im fluxes. This resulted in a 
list of 39 objects for which the photometric data are given 
in Appendix B. The characteristics of these 39 objects are 
given in Table 1. An additional remarkable RSG, NML Cyg, 
will be examined in the next subsection. 

In order to compare the mass-loss rate of these RSGs 
with the de Jager prescription, we estimate M from the 
infrared excess at 60 (im with the formulation of Jura and 
collaborators (JK90): 

M 60AJm = 1.7 x 10- 7 V 15 Dl pc f 60 (A 10 /£ 4 ) 1/2 , 

where M 6 o^ m is in M Q yr _1 , V15 is the wind expansion 
velocity normalized to 15 kms -1 , D^ pc is the distance in 
kpc, /go is the IRAS flux at 60 (im in Jy (assumed to domi- 
nate the photospheric flux at this wavelength), Am is equal 
to the mean wavelength A m of the object spectral energy 
distribution, normalized to 10 /im, and L4 is the luminosity 
in units of 10 4 L Q . This formula is valid for a gas-to-dust 
mass ratio of 200, which corresponds to complete conden- 
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Table 1. Properties of the RSGs. The first 27 stars are from LM05, while the following 12 stars are from JBOO. The column 
marked V is the assumed wind speed. The column marked A is the flux averaged wavelength. L^ ev is the luminosity given 
in LM05 (from their Mboi) and L p hot is the luminosity obtained by integrating the UBVIJHKL + IRAS photometry. M 
is derived with Jura's formula, where the used luminosity is L p hot- 
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rn — M 




f 1 2 / f K 


V 


D 


,/*60 






-^phot 




M 




(mag) 


(mag) 




(km s ) 


(kpc) 


(Jy) 


H 


(L©) 


(L©) 


(M yr- 1 ) 


V589 Cas 


11.50 


0.78 


0.156 


14 


2.00 


3.61 


1.72 


52 000 


35 000 


5.0 


x 10~ 7 


BU Per 


11.40 


0.66 


0.435 


14 


1.90 


5.23 


2.20 


58 000 


38 000 


7.4 


x 10~ 7 


SU Per 


11.40 


0.66 


0.241 


19 


1.90 


6.87 


1.77 


90 000 


85 000 


7.7 


x 10~ 7 


RS Per 


11.90 


0.56 


0.417 


20 


2.40 


9.93 


2.28 


144 000 


95 000 


2.0 


x 10 -6 


S Per 


11.39 


0.66 


1.226 


20 


1.90 


40.59 


3.67 


81000 


86 000 


6.8 


x 10~ 6 


V441 Per 


11.40 


0.66 


0.156 


16 


1.90 


3.54 


1.69 


66 000 


50 000 


4.2 


x 10~ 7 


YZ Per 


11.40 


0.66 


0.291 


16 


1.90 


5.28 


1.86 


48 000 


55 000 


6.5 


x 10~ 7 


W Per 


11.40 


0.66 


0.495 


16 


1.90 


14.87 


2.55 


54 000 


56 000 


2.1 


x 10~ 6 


BD+57647 


11.40 


0.66 


0.351 


14 


1.90 


6.47 


2.24 


80 000 


37000 


9.3 


x 10~ 7 


NO Aur 


10.70 


0.47 


0.146 


18 


1.38 


5.12 


1.62 


67000 


73 000 


2.9 


x 10~ 7 


a Ori 


5.57 


0.18 


0.177 


15 


0.130 


299.00 


1.64 




56 000 


1.5 


x 10~ 7 


TV Gem 


10.70 


0.66 


0.294 


19 


1.38 


6.06 


1.69 


100 000 


84000 


3.5 


x 10~ 7 


BU Gem 


10.70 


0.66 


0.248 


19 


1.38 


10.50 


1.62 


83 000 


86 000 


5.9 


x 10~ 7 


V384 Pup 


13.00 


0.57 


0.260 


18 


4.00 


2.76 


1.84 


37000 


75 000 


1.4 


x 10~ 6 


CK Car 


11.70 


0.55 


0.531 


22 


2.20 


13.98 


2.11 


161 000 


123 000 


2.1 


x 10~ 6 


V602 Car 


11.60 


0.48 


0.596 


22 


2.10 


12.40 


2.61 


105 000 


124 000 


1.9 


x 10~ 6 


V396 Cen 


11.60 


0.75 


0.201 


23 


2.10 


4.98 


1.73 


164 000 


140 000 


6.2 


x 10~ 7 


KW Sgr 


12.40 


0.92 


1.218 


27 


3.00 


18.39 


2.81 


363 000 


228 000 


5.6 


x 10~ 6 


NR Vul 


11.80 


0.94 


0.590 


21 


2.30 


12.28 


2.26 


224 000 


111000 


2.2 


x 10~ 6 


BI Cyg 


11.00 


0.93 


0.671 


22 


1.58 


51.23 


2.67 


226 000 


123 000 


4.6 


x 10~ 6 


KY Cyg 


11.00 


0.93 


0.702 


22 


1.58 


50.74 


3.15 


272 000 


138 000 


4.9 


x 10~ 6 


RW Cyg 


10.60 


1.22 


0.481 


23 


1.32 


60.69 


1.96 


144 000 


145 000 


3.2 


x 10~ 6 


/i Cep 


9.70 


0.69 


0.361 


20 


0.87 


127.00 


1.69 


340 000 


410 000 


1.4 


x 10~ 6 


V354 Cep 


12.20 


0.63 


0.566 


18 


2.75 


8.01 


2.89 


369 000 


76 000 


2.4 


x 10~ 6 


V355 Cep 


12.20 


0.63 


0.292 


14 


2.75 


3.27 


2.38 


94 000 


37000 


1.0 


x 10~ 6 


PZ Cas 


11.90 


0.69 


1.217 


30 


2.40 


96.48 


3.57 


212 000 


193 000 


2.6 


x 10~ 5 


TZ Cas 


11.90 


0.69 


0.541 


19 


2.40 


9.47 


2.38 


98 000 


83 000 


2.0 


x 10~ 6 


EV Car 


13.13 


0.51 


0.745 


39 


4.20 


25.87 


2.57 




675 000 


1.3 


x 10 -6 


HS Cas 


12.00 


0.83 


0.275 


17 


2.50 


3.51 


1.92 




59 000 


7.5 


x 10~ 7 


XX Per 


11.14 


0.66 


0.492 


16 


1.69 


4.23 


2.03 




50 000 


4.3 


x 10~ 7 


KK Per 


11.14 


0.66 


0.117 


16 


1.69 


2.23 


1.59 




50 000 


2.0 


x 10~ 7 


AD Per 


11.14 


0.66 


0.164 


15 


1.69 


2.85 


1.66 




42 000 


2.7 


x 10~ 7 


PR Per 


11.14 


0.66 


0.131 


14 


1.69 


2.37 


1.55 




34000 


2.3 


x 10~ 7 


GP Cas 


11.40 


0.66 


0.200 


15 


1.90 


4.45 


2.01 




43 000 


5.9 


x 10~ 7 


VY CMa 


10.28 


0.47 


6.959 


47 


1.14 


1453.00 


7.77 




295 000 


1.6 


x 10~ 4 


a Sco 


6.34 


0.10 


0.165 


17 


0.185 


115.50 


1.73 




71000 


1.2 


x 10~ 7 


VX Sgr 


10.98 


0.52 


2.334 


25 


1.57 


262.70 


4.14 




343 000 


2.0 


x 10 _s 


Case 49 


11.70 


0.87 


0.155 


15 


2.19 


6.58 


2.01 




42 000 


1.1 


x 10~ 6 


U Lac 


12.20 


0.63 


0.685 


23 


2.75 


9.04 


2.51 




147 000 


2.3 


x 10~ 6 



sation of refractory elements for solar abundances, and we 
adopt this ratio here. 

The wind speed in Jura's formula is known for only 
a few RSGs (see Table 2). For all others, we choose to 
scale the value for this speed with the luminosity, and we 
adopt V oc L 35 following predictions of dust-driven wind 
modelling (Habing et al. 1994). For the RSGs of the solar 
neighbourhood, we find that V ~ 20 x (i/10 5 ) 35 km s" 1 
roughly represents the observed values. For the LMC RSGs 
considered below, we find V ~ 14 x (L/10 5 ) - 35 km s _1 (see 
Appendix C for details). 

Concerning the luminosity and mean wavelength, these 
quantities are calculated by integrating over the spectral en- 



ergy distribution. We use the UBVIJHKL magnitudes and 
IRAS fluxes and correct the magnitudes for interstellar red- 
dening following the extinction law of Cardelli et al. (1989). 
Between the central wavelengths of these spectral bands, in- 
terpolation is done linearly on the relation of log(A f\) as 
a function of log(A). The accuracy of the derived luminosi- 
ties (listed as £ p hot in Table 1) depends on several factors. 
Once the distance is fixed with the distance modulus, fac- 
tors like photometric uncertainties, variability or errors on 
Eb-v introduce uncertainty. For example, adopting JHK 
magnitudes smaller by 0.2 mag results in larger luminosity 
by 11%, and decrease of M by & 7%. Similarly, increas- 
ing Eb-v by 0.25 mag implies an increase in L by ^ 40%, 
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and a decrease in M by ^ 16%. These variations are smaller 
than the accuracy of the Jura's formula. 

In Table 1, we have included the luminosities from LM05 
(column Llcv for the 26 objects originating from their 
study). Many determinations of luminosity are in very good 
agreement between theirs and ours. However, there are also 
a few significant discrepancies, up to a factor of 5 for V354 
Cep. Because we use the same distance and reddening, and 
almost the same photometry, the discrepancy must origi- 
nate from the method for deriving L, but we could not trace 
exactly where the differences come from. Consequently, we 
used our luminosities for plotting M versus L (Fig. 2). 

It can be seen in Fig. 2 that, if we omit four exceptions 
labelled 1,2,3 and 7, the majority of points are clustered 
near the de Jager prescription line. For \ogL < 5.2, the de 
Jager line slightly overestimates the mass-loss rates (by a 
factor ~ 1.5), while for \ogL > 5.2, the de Jager line passes 
below the stars PZ Cas and VY CMa, but above KW Sgr 
and VX Sgr, and increasing the dJ88 rate by 1.5 would not 
improve the situation. 

The stars a Ori, a Sco and \x Cep (labelled 1, 2, 3 re- 
spectively) present a low M here, but observations of their 
circumstellar gas suggest that they are plausibly depleted 
in grains (see next section). We note that our estimate of M 
for \i Cep, based on the 60 /xm flux, is in good agreement 
with that of de Wit et al. (2008) derived from imaging at 
24 /jm, after appropriate scaling of the outflow velocity. 

There is also one star at very high luminosity: EV 
Car (labelled 7), with L = 6.7 x 10 5 L Q and M =1.3 x 
10 -5 M Q yr _1 . The huge luminosity is due to the distance 
modulus given by Humphreys (1978), (m — M) = 13.3, cor- 
responding to 4.2 kpc. This puts the star far beyond the Car 
OBI association. If we had taken the distance of this asso- 
ciation (2.5 kpc, adopted by JK90), we would have found 
L = 2.4 x 10 5 L and M = 4.6 x 10~ 6 M Q yr" 1 , and the 
location of the star in Fig. 2 would be very close to the de 
Jager line. But to our knowledge, there is no observational 
evidence for EV Car to be member of Car OBI. 

3.2. Circumstellar gas of RSGs in the solar neighbourhood 

Only a few RSG winds have been detected through observa- 
tion of the circumstellar gas. Methods of detection include 
measurements of the Hi 21cm emission (details on observa- 
tions and mass-loss determination can be found in Bowers 
& Knapp 1987), millimetric emission of CO (Olofsson 2004; 
Ramstedt et al. 2008), submillimctcr emission of C I (Kcenc 
et al. 1993; Huggins et al. 1994), emission of scattered 
resonance lines of Nal and Ki (Mauron 1990; Guilain & 
Mauron 1996; Gustafsson et al. 1997), and observation of 
the Hn region around Antares (Kudritzki & Rcimcrs 1978; 
Reimers et al. 2008). 

For eight stars, we have collected in Table 2 the mass- 
loss rates obtained with these methods, giving detailed ref- 
erence for adopted parameters (D and Eb-v) use d to de- 
rive luminosities. When M was obtained in original papers 
with a different distance, we assumed M ocD 2 . The expan- 
sion velocities V are from Huggins et al. (1994) for a Ori, 
Woodhams (1993) for S Per and PZ Cas, and from Kemper 
et al. (2003) for VY CMa, VX Sgr, NML Cyg. Finally, for 
\i Cep, the situation is unclear (for details see Bernat 1981, 
the discussion in Le Borgne & Mauron 1989, and Kemper 
et al. 2003) and we have adopted 20 kms -1 . 




4 4.5 5 5.5 

log (Luminosity) 



Fig. 2. Comparison of de Jager rate (plotted for T e s = 3500 
K and 4000 K) with mass-loss rates of 39 RSGs listed in 
Table 1. The mass-loss rates are based on the study of dust 
and obtained with Jura's formula. The labels indicate a few 
wellknown RSGs as follows: 1: a Ori; 2: a Sco, 3: /i Cep, 4: 
S Per, 5: PZ Cas, 6: VY CMa, 7: EV Car (at 4.2 kpc); 8: VX 
Sgr. The small cross represents EV Car if a smaller distance 
of 2.5 kpc is adopted. At the upper left corner, error bars 
represent uncertainties of ± 30 percent in luminosity, and 
a factor of ± 2 in M. 



Several remarks have to be made on individual objects. 
For the most observed star (a Ori), there is good agreement 
among the rates based on H I, C I and K I (Bowers & Knapp 
1987; Huggins et al. 1994; Glassgold & Huggins 1986; Plez 
& Lambert 2002). When the distance given by Hipparcos 
is adopted (parallax 7.63 ± 1.64 mas, giving 130 pc ± 30 pc), 
one obtains M « 1 to 2 x 10~ 6 M Q yr _1 . The distance pro- 
posed by Harper et al. (2008) is 197 ± 45 pc and would im- 
ply M ~ 3x 10 -6 M Q yr _1 . For this star, the interstellar ex- 
tinction Ay is not very well determined, being between 0.3 
and 0.8 mag according to Lambert et al. (1984). We have 
chosen the average of these values implying L — 56000 L , 
but the main uncertainty comes from the error on parallax, 
and L is known within ± 50 percent. 

In the case of VX Sgr, the distance has recently been 
measured by Chen et al. (2007) with SiO maser proper mo- 
tions (D = 1.57 ± 0.27 kpc), and this distance means that 
the star belongs to the Sgr OBI association. The average 
interstellar extinction to the OB stars of this association 
is Eb~v = 0.5 (from data of Humphreys 1978), and the 
luminosity is 3.4 x 10 5 L Q . 

In the case of VY CMa, a distance of 1.14 ±0.10 kpc 
has been obtained by Choi et al. (2008) through astro- 
metric monitoring of H2O masers. This distance is smaller 
than often used. For example, Dccin et al. (2006) have 
adopted 1.50 kpc for their detailed study of circumstel- 
lar CO lines. The ejection rate obtained by Decin et al. is 
highly variable. The episode with the highest mass-loss rate 
(~ 3 x 10~ 4 M yr _1 ) has a duration of only 100 yr. 
The average mass-loss rate is about 7 x 10~ 5 M Q yr~ 1 . 
Scaled to the distance of 1.14 kpc, the average rate becomes 
4x 10~ 5 Moyr -1 . The new smaller distance means that 
the star is less luminous than previously believed. We find 
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NML Cyg 




4 4.5 5 5.5 

log (Luminosity) 

Fig. 3. Mass-loss rates (from observations of the gas) 
versus luminosity for the eight RSGs of Table 2. The lines 
are the de Jager rates for 4000K (solid line) and 3500K 
(dotted line). 



mass-loss rates of red supergiants 




-7 -6 -5 -4 -3 

log(Mass loss rate from dust) 



Fig. 4. Mass-loss rates from observation of the gas versus 
mass-loss rate from observation of the dust (the 60 ^m emis- 
sion) for the eight RSGs of Table 2. 



L = 2.9 x 10 5 L©. Because the star is extremely red (the flux 
averaged wavelength is around 8/im, compared to ~ 2um 
for the majority of objects), this luminosity does not de- 
pend much on the choice of the extinction: for -Eb-v m 
the range to 1, the luminosity is between L = 2.7 and 
3.5 x 1O 5 L . 

The points representing these eight stars with measured 
gas mass-loss rates as a function of luminosity are shown in 
Fig. 3. The error bars correspond to ±30 percent in lumi- 
nosity and a factor of ± 2 of uncertainty on the mass-loss 
rate. The actual uncertainties on M are probably larger 
than that but difficult to quantify. The points are not very 
far from the de Jager line. Betelgeuse and Antares agree for 
M w 10~ 6 Moyr" 1 at a luminosity of - 65000 L . SPer 
is a factor of 4 above the line and \i Cep is a factor of 10 
below the line. It can be seen that on average the de Jager 
prescription agrees reasonably well with measurements of 
the gas mass-loss rates. 

Since Meoum was used in the previous section, we can 
compare it, for these stars, with the mass-loss rates derived 
from observation of the gas M gas - This is done in Fig. 4 and 
Table 2. Figure 4 shows that the agreement between the 
two rates is satisfactory for some stars (S Per, NML Cyg, 
VX Sgr), but there are also differences of a factor of ~ 4 
for VY CMa and u Cep. These differences may arise from 
the fact that different regions of the envelope are probed 
when comparing different techniques. For a Ori and a Sco, 
M e0finl is too low by a factor of <~10, plausibly due to in- 
complete dust formation. So, there is a trend in Fig. 4 that 
for low M, grain condensation may be incomplete, but ob- 
viously many additional sources should be examined with 
independent estimates of the dust mass-loss rate and the 
gas mass-loss rate to confirm this trend. 

3.3. Infrared measurements in the Large Magellanic Cloud 

The mass-loss rate of RSGs in the Large Magellanic Cloud 
(LMC) has been investigated in several works, with vari- 



ous data and various methods. All works are based on the 
circumstellar dust emission. Reid et al. (1990) used IRAS 
fluxes and the Jura's formula, van Loon et al. (2005) ex- 
ploited MSX, ISO and IRAS fluxes (their Table 3) and anal- 
ysed the spectral energy distribution with the dust radiative 
transfer code DUSTY, which provides M. More recently, 
Groenewegen et al. (2009) exploited Spitzer results in ad- 
dition to all previously obtained data. In particular, Spitzer 
data include photometry from 3.6 to 24 um and spec- 
troscopy from 5.2 to 37 /im. Groenewegen et al. (2009) fit 
the data also with dust radiative transfer models. (The work 
of Bonanos et al. (2010) also treats mass loss of RSGs in 
the Magellanic Clouds and is considered in the Discussion.) 

The mass-loss rates given in the works cited above can- 
not be compared directly, because their assumed param- 
eters, such as the expansion velocity of the wind or the 
gas-to-dust mass ratio are different. We assume here that 
the LMC RSGs have a metallicity lower than that of the 
Sun so that complete condensation of refractory metals 
leads to a gas-to-dust mass ratio tp — 500. We also as- 
sume that all RSGs have a wind velocity v(L) depend- 
ing on the luminosity L as found in Appendix C i.e. 
v(L) = 14(L/10 5 ) 35 kms -1 . Consequently, we have multi- 
plied the M of Reid et al. (1990) by 2.5 x v(L)/15 (because 
the Jura's formula used by Reid et al. supposes = 200 
and a wind speed of 15 kms -1 ). We have multiplied the 
M of van Loon et al. (2005) by v(L)/v' where v' is their 
published wind velocity provided by DUSTY. Finally, we 
have multiplied by 2.5 x v(L)/10 the M of Groenewegen 
et al. (2009) because they adopt ip — 200 and a wind speed 
of 10 kms -1 . 

The results are in Table 3, 4 and 5. As far as possible, we 
did not include AGB stars in these tables. More precisely, 
we followed Groenewegen et al. (2009) by considering that 
all objects with I ^ 1.2 x 10 5 L© are genuine RSGs, and 
that objects with L between ~ 5 x 10 4 and 1.2 x 10 5 L Q can 
be considered as RSGs only if their amplitude of variation 
in the optical / band was less than 0.45 mag. This selection 
implied that 7 objects in the Groenewegen sample were kept 
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Table 2. Characteristics of RSGs for which a mass-loss rate has been measured by observing the circumstellar gas. 



Name 


D 


-£/B— V 


V 


■Feo 


Am 


iphot 


M QQfj,m 


Mgas 


Method 


References 




fkpcl 


[mael 

L oj 


[km s _1 ] 


Jv 


\uxn\ 

If-*' \ 




[10~ 6 M 


p, vr _1 l 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


a Ori 


0.130 


0.18 


15 


299 


1.64 


56 000 


0.15 


1.5 


Hi Ci Ki 


1 


2 


3 


a Sco 


0.185 


0.10 


17 


115 


1.73 


71000 


0.12 


1.0 


Hiireg. 


1 


4 


5 


S Per 


1.90 


0.66 


20 


41 


3.67 


86 000 


6.8 


7.5 


CO 


6 


6 


7 


PZ Cas 


2.40 


0.69 


30 


96 


3.57 


193 000 


26 


8.3 


CO 


8 


8 


7 


VY CMa 


1.14 


0.47 


47 


1453 


7.77 


295 000 


160 


40 


CO 


9 


10 


11 


NML Cyg 


1.74 


1.20 


34 


1030 


9.80 


320 000 


210 


140 


CO 


12 


13 


14 


VX Sgr 


1.57 


0.52 


25 


263 


4.14 


340 000 


20 


11 


CO 


15 


16 


17 


fj, Cep 


0.87 


0.69 


20 


127 


1.69 


410 000 


1.4 


5.0 


Ki Nai 


6 


6 


18 



Notes: Col. (11) = reference for distance. Col. (12) = reference for _Eb-v- Col. (13) = reference for mass-loss rate M gas . The IRAS 
data of NML Cyg were taken from Schuster et al. (2009) . 

References: (1) Hipparcos parallax; (2) Lambert et al. 1984; (3) Huggins et al. 1994; (4) assumed; (5) Reimers et al. 2008; (6) 
LM05; (7) Woodhams 1993; (8) JB00; (9) Choi et al. 2008; (10) assuming E B -v = 0.47 per kpc; (11) Decin et al. 2006 & see text; 
(12) Massey & Thomson 1991; (13) Blocker et al. 2001; (14) Knapp et al. 1982; (15) Chen et al. 2007; (16) Humphreys 1978; (17) 
Knapp et al. 1989; (18) Guilain & Mauron 1996. 



(marked with "s" after their names in the first column of 
Table 5) and 9 were rejected. 

Table 3. Red supergiants of Reid et al. (1990). The 
columns provide the object name, J2000 coordinates as for 
2MASS names, luminosities in L , mass- loss rates in units 
of 10 -6 M Q yr~ 1 , and the ratios r with mass-loss rates 
from Groenewegen et al. (2000) for common objects. 



Name 


a, 8 (J2000) 


L 


M 


r 


HV 12501 


045541.83- 


-692624.3 


132 000 


157 




HV 12420 


045731.52- 


-700859.2 


110 000 


79 




HV 2255 


045743.32- 


-700850.3 


251 000 


196 


37 


HV 888 


050414.14- 


-671614.4 


328 000 


134 


37 


HV 894 


050533.50- 


-703347.0 


106 000 


86 




HV 2360 


051246.37- 


-671937.9 


161 000 


135 




HV 916 


051449.74- 


-672719.8 


148 000 


96 


37 


HV 2532 


052627.39- 


-691055.9 


113 000 


49 




HV 963 


052734.35- 


-665330.0 


100 000 


56 


31 


HV 5854 


052815.41- 


-665852.7 


97 000 


23 




HV 2586 


052934.43- 


-665528.3 


71 000 


50 




HV 996 


053235.61- 


-675508.9 


127 000 


172 


26 


HV 12437 


053307.61- 


-664805.7 


82 000 


65 




HV 2700 


053518.92- 


-670219.6 


89 000 


42 




HV 1004 


053625.46- 


-665538.4 


141 000 


110 





Tables 3, 4 and 5 include names, J2000 coordinates, 
luminosities, mass-loss rates (scaled as explained above and 
in units of 10 _6 M Q yr _1 ). The last column in Table 3 and 
4 are the ratio of the M from Reid et al. (or from van 
Loon et al. ) and the M from Groenewegen et al. (2009) 
for objects belonging to both samples. 

Figure 5 shows these scaled M as a function of L from 
the three works mentioned above. We have also shown the 
de Jager rate for 4000 and 3500 K exactly as for previous 
plots. It can be seen immediately that there is a difference 
between the recent rates from Groenewegen sample on one 
hand (the median rate is 2.3 x 10~ 6 Mgyr" 1 , 30 RSGs 
represented by circles), and, on the other hand, the Reid 



Table 4. Red supergiants of van Loon et al. (2005). Same 
columns as in Table 3. 



Name a, S (J2000) L M 



WOH-G64 


045510.48- 


-682029.8 


490 000 


1200 


1.6 


HV 12501 


045554.83- 


-692624.3 


145 000 


15 




HV 888 


050414.14- 


-671614.4 


269 000 


37 


10 


HV 2360 


051246.37- 


-671937.9 


129 000 


30 




HV 916 


051449.74- 


-672719.8 


148 000 


35 


13 


HV 2561 


052828.86- 


-680707.9 


179 000 


49 


16 


HV 5870 


052903.48- 


690646.2 


95 000 


23 


8 


SP77-46-44 


052942.21- 


-685717.4 


295 000 


52 


37 


HV 986 


053109.28- 


-672554.9 


132 000 


20 




HV 996 


053235.61- 


-675508.9 


117 000 


63 


10 



sample (median 9.0 x 10~ 5 M Q yr _1 ; 15 crosses) or the 
van Loon sample (median 3.6 x 10~ 5 Moyr -1 ; 10 small 
squares). One can also consider the objects in common (see 
Table 3 and 4): these objects indicate that the Reid mass- 
loss rates are ~ 35 times larger than those of Groenewegen, 
and that the van Loon mass-loss rates are ~ 10 times larger 
than those of Groenewegen. One notes however that, for 
WOH-G64, the M from van Loon and Groenewegen agree 
reasonably well (within a factor 1.6). In contrast, for SP77- 
46-44 (alias W60 A27), they disagree by a factor 37. 

There are two assumptions in Reid et al. (1990) which 
can explain (at least in part) why their rates are so high. 
First, there is probably an overestimate of the flux averaged 
wavelength: Reid et al. (1990) take Aio = 1, which is not 
reached even by NML Cyg (see Table 2). For the optical 
RSGs making the Reid sample, Aio is likely close to 0.2 
only. Secondly, Reid et al. (1990) overestimate the f§o/ /2s 
IRAS flux ratio used to obtain f^. They take 0.50 instead 
of ~ 0.20. A ratio of 0.23 ±0.03 is reached only by dusty 
Galactic RSGs (VX Sgr: 0.19, VY CMa: 0.22, PZ Cas: 0.24, 
NML Cyg: 0.26). So, the Reid et al. mass- loss rates have to 
be decreased by a factor ^12. 

Concerning the van Loon rates, the reasons why they are 
generally larger than those of Groenewegen are not clear. 
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Table 5. Red supergiants of Groenewegen et al. (2009), 
with names, J2000 coordinates, luminosities and mass-loss 
rates in units of 10 -6 M yr _1 . 



Name 


a, 5 (J2000) 


L 


M 


HV 99^8 

11 V L^OU 


n a a noo a ft 
U44yzz.4o — 


-oyz4o4.4 


1 9n nnn 

1ZU uuu 


9 8 


MSX IMC 1 189 

IVlkJW UlUVj 110^7 


n A K £n*3 H7 


oyzyiz. ( 


1 3Q 000 

lot/ uuu 


6.7 


WOH Cif\d 


U4oolU.4o- 


DozUzy.o 




7R4 


MSX T,MC 1904 

IVlkJW UlUVj liiU4 


n /i £ £ i ft. 


ftm ni o 1 

oyiyiz.i 


904 000 

Zj\J l ± uuu 


0.81 


msx T.Mn ms c 

1V1QYV J-JlvlV^ lOlO o 


U4oooo.o4- 


ftno icri o 


fil 000 
Ul uuu 


1 

l.O 


HV 99 c i c i 

11 V 


U4o ( 4o.ol- 


- / UUooU.o 


1 4*3 000 
l^O uuu 


fJ.O 


HV 888 
11 V ooo 


C\KC\A 1 A 1 Q 

UoU414. lo- 


ft71 CI A Q 

/ 1014.O 


QQ7 000 
OO I uuu 


o.u 


HV 916 


ne; 1 /I /I O TO 

Uol44y. / z — 


otZt iy. / 


9^0 000 

ZiOU uuu 


2.6 


WDH S9(S4 

VV Wll Oiul 


Uoz4iy.oz — 


-oyoo4y.o 


91 ^ 000 

ZltJ uuu 


4 ^ 


MSX LMC 540 


ne:ofti i 
Uozull.oO- 


ftft 1 Oil 1 

DDlzll.l 


1 S9 000 

lOZ/ uuu 


2.0 


msx imp 

IVlQYv 1j1V1\_> O I O 


Uozuzz.lo— 


DDzlzo.o 


1 fi^ 000 
1UO uuu 


n 9n 


MSX T,MC 589 


AKOftQ/1 on 

U0zuo4.oU- 


fioci /in n 


9fiQ 000 

£i\j<j uuu 


1.7 


HV Qfi^ « 
ii v yuo o 


Uoz f o4.oo- 


DDOOOU.U 


1 1 ^ 000 

1 10 uuu 


1 8 

l.O 


HV 9561 

11 V 4j<J\J 1 


nC^OQOQ Cft 


DoU ( U / .0 


1 000 

lot/ uuu 


3.0 


HV ^870 « 

11 V OO 1 U o 


UozyUo.4c>- 


ftnnft /ic o 
-DyUD4o.Z 


7Q 000 

/ a UUU 


9 8 


Wfin A 97 


nenn /I o Ol 

Uozy4z.zl — 


ft Q 7 Q 
OOO i ii .O 


OQC] noo 

ZiOU uuu 


1.4 


1V1 -\ 1j 1V1 OIU 


A^QflOn ft7 


DDOoUl.o 


1 Q7 000 
iy i uuu 


1.1 


MSX f,MC 587 


UoolU4.1o- 




91 R 000 

ZilO uuu 


2.3 


WfiO F)99 


nroi in ft ,1 

053110.64- 


ft ft O C O 1 ft 

-663531.6 


1 fi^ 000 
1UO uuu 


u.uu 


yqy t yp qqq 


053136.81- 


-663007.6 


9^^ 000 
zoo uuu 


9 9 


HV 996 


053235.61- 


-675508.9 


132 000 


6.6 


MSX LMC 791 s 


053524.52- 


-690403.4 


84 000 


0.89 


MSX LMC 870 


053528.32- 


-665602.4 


218 000 


2.3 


MSX LMC 891 


053555.22- 


-690959.4 


265 000 


3.1 


W60 A72 s 


053833.97- 


-692031.7 


115 000 


84 


W60 A78 


053932.33- 


-693450.0 


223 000 


0.79 


MSX LMC 897 


054043.75- 


-692158.1 


245 000 


1.0 


MSX MC 939 s 


054048.50- 


-693336.1 


109 000 


4.4 


HV 1017 s 


054059.20- 


-691836.2 


102 000 


4.3 


HV 2834 


054413.73- 


-661644.5 


215 000 


1.1 



With Spitzcr data, small infrared excesses become measur- 
able and are more accurate, while this was not the case 
with MSX, ISO or IRAS data available for the van Loon et 
al. work. Therefore, there may be a bias in favour of large 
M in the van Loon sample compared to the Groenewegen 
sample. Another possibility is that the higher spatial reso- 
lution of Spitzer allows to avoid source confusion and con- 
taminating sources. This may be in particular the case for 
SP77-46-44. For this source, the IRAS fluxes at 12 and 25 
jum are 0.26 and 0.18 Jy, while Spitzcr photometry (from 
Bonanos et al. 2009) indicates 0.19 Jy at 12 /xm (through 
interpolation from the 8 (im and 24 fluxes) and 0.11 Jy 
at 25 yum. So these fainter Spitzer fluxes suggest a smaller 
infrared excess and may go some way to explain the large 
discrepancy for this object. 

Finally, if we consider the rates by Groenewegen et 
al. (2009) and their positions with respect to the de Jager 
lines (see Fig. 5), one can note a rather good agree- 
ment for 12 objects with luminosities lower than 1.6 x 
10 5 L (with the exception of W60-A72 at L = 115000 L Q , 
M = 8.4 x 10~ 5 M yr _1 ; in Fig. 5, this object is the circle 
located in the middle of Reid's plus signs). But at larger lu- 
minosities, the Groenewegen rates are below the dJ88 lines 
by an average factor ~ 8 (13 objects), with the exception 
of WOH-G64 where the de Jager rate is ~ 10 times lower 
than observed. It is also worthy to note that at high lumi- 
nosities, two objects from the van Loon sample are quite 
close to the dJ88 line. 
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Fig. 5. Mass-loss rates plotted as a function of luminosity 
for the LMC RSGs. Three samples are plotted: Reid et 
al. (1990) (15 small crosses), van Loon et al. (2005) (10 
small squares), and Groenewegen et al. (2009) (30 circles). 
The de Jager rates are represented as in previous figures by 
solid and dotted lines for 4000 K and 3500 K, respectively. 



4. Discussion 

The goal of this paper is to see whether the de Jager pre- 
scription was still applicable in view of new results obtained 
since 1988 and concerning the RSGs and their winds. One 
important fact is that for only eight RSGs in the solar 
neighbourhood is the direct detection of the circumstcllar 
gas available. In all other cases, only dust can be seen. We 
find that for these 8 RSGs, the mass-loss rate estimated 
from observations of the gas is in most cases within a factor 
of 4 of the de Jager prescription for L < 2 x 10 5 L (Fig. 2). 
For higher luminosity, it is a factor of 10. These devia- 
tions have to be considered in light of the following char- 
acteristics of the RSG winds. The winds are changing with 
time (timescale ~ 10 2 - 10 4 years), and the rate of mass 
loss and the expansion velocity vary. This is proved by 
high resolution spectroscopy of the circumstellar absorption 
lines which display discrete velocity components (on \x Cep, 
Bernat 1981), by direct imagery (a Ori, Plez and Lambert 
2002) or by study of the CO millimetric lines (VY CMa, 
Decin et al. 2006). Also, the winds are often not spherically 
symmetrical, as revealed by some CO line profiles (PZ Cas, 
Woodhams 1993) or direct imagery (HST imaging of NML 
Cyg, Schuster et al. 2006). Given these properties, the mea- 
surements of M with often limited data and simple mod- 
elling are necessarily poor in accuracy and do not reflect an 
average mass-loss rate. Therefore, the agreement with the 
de Jager prescription seems acceptable. 

For other, much more numerous RSGs, the gas is not 
detectable and only dust can be seen through its infrared 
radiation. When the Jura's formula is adopted for deriving 
M, the resulting M are also in agreement with de Jager 
lines, except for some rare stars like \x Cep which has a 
very low dust content. In general, the agreement is within 
a factor of 4 (Fig. 2). Part of this deviation may be due 
to a natural scatter of M at a given luminosity (see the 
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characteristics of winds mentioned above) and another part 
may be due to the use of the Jura's formula. 

4.1. Verification of the Jura's formula 

In order to check the Jura's formula, we have considered 
the work by Verhoelst et al. (2009) who exploit ISO-SWS 
spectra (in addition to the full spectral energy distribu- 
tion) and use state-of-the-art dust modelling. In particular, 
they include several types of grains such as melilite, olivine, 
alumine, carbon and FeO. Together with results concern- 
ing the dust condensation sequence in RSGs, they derive 
mass-loss rates. We took into account their assumptions 
on distance, expansion velocity and gas-to-dust mass ratio 
and we applied appropriate scaling, so that we can compare 
their results to ours (Meo^m in Table 2). Wc find that their 
M are 0.57, 0.88, 1.20 and 1.9 times greater than ours, for 
a Ori, S Per, PZ Cas and /x Cep, respectively. In conclu- 
sion, the agreement with Jura's formula is at most a factor 
of 2, which is very good, and there is no particular reason 
to have doubts on this formula. 

4.2. Comparison with other mass-loss prescriptions 

In Sect. 2, we mentioned in Fig. 1 a few mass-loss prescrip- 
tions. After our analysis of previous sections, we can recon- 
sider this point. Concerning the RSGs in the solar vicinity, 
all observed mass-loss rates (from dust or gas observations) 
are well below the rate by Salasnich et al. (1999). For exam- 
ple, a Ori and a Sco are located at L rts 60 000 L©, and their 
M are ~ 1.5 x 10~ 6 M© yr _1 . The Salasnich rate would give 
3.4 x 10~ 5 M© yr _1 , 20 times larger. Because they are based 
on mass- loss rates of Reid et al. which we have shown are 
overestimated, the Salasnich rate and the Feast rate are 
too large. The Reimers rate is as good as the dJ88 rate for 
objects with L ^ 2 x 10 5 L©, but this is not the case for 
higher luminosity objects like NML Cyg or VY CMa. 

4.3. Comparison with the van Loon formalism 

It is appropriate to consider here the mass-loss formu- 
lation of van Loon et al. (2005). These authors derived 
a prescription for mass-loss rates of AGB stars and red 
supergiants. It is based on the study of evolved stars in 
the LMC, but is also valid for Galactic AGB stars and 
RSGs. An important feature of this prescription is that it 
is applicable only to dusty stars. Unfortunately, the dusty 
character of a circumstellar envelope is not possible to 
be derived from the basic stellar parameters, i.e. L, M 
and T c g. Nevertheless, one can compare the results of this 
prescription to the de Jager rate, van Loon et al. (2005) 
derived effective temperatures from spectral classifications. 
Luminosities and mass-loss rates were determined by 
modelling the spectral energy distribution. The van Loon 
prescription is: 

M = 2.50 x 10~ 5 (L/10 5 ) 1 - 05 (T/3500)- 6 - 3 . 

This prescription is shown in Fig. 6 along with the de 
Jager lines. Both are drawn for T e s of 4000, 3500 and 
3000 K. van Loon et al. (2005) have found that their recipe 
overestimates mass-loss rates for Galactic "optical" RSGs. 
We can see that clearly in Fig. 6, where the van Loon lines 
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Fig. 6. Mass-loss rate prescription of van Loon et al. (2005) 
and de Jager (1988) with M plotted as a function of lumi- 
nosity and for three values of T c g. The de Jager lines are 
solid, dotted and dashed for 4000, 3500 and 3000 K respec- 
tively. 

are well above the de Jager lines on which we found that 
most of the Galactic RSG are located (Fig. 2 and Fig. 3). 
However, the de Jager rate and the van Loon rate fairly 
agree with each other for the highest luminosities. Also, it 
is worth noting that the van Loon et al. prescription offers 
a significant spread of M for a given luminosity. The spread 
allowed by the dcJ88 prescription is clearly very small. 

We have collected in Table 6 the mass-loss rates for 
the eight stars possessing an M gas based on observation of 
the gas (from Table 2). In this Table, the T e g values are 
taken from Lambert et al. (1984) for a Ori, from LM05 for 

5 Per, PZ Cas, and /z Cep, from van Loon et al. (2005) 
for VY CMa, NML Cyg and VX Sgr, and from Kudritzki 

6 Reimers (1978) for a Sco. Column 5 of this Table gives 
the de Jager rates, and column 6 gives the van Loon rates 
(in column 6, the numbers in parenthesis concern "opti- 
cal" RSGs and are just given for information). The ratios 
between van Loon rates and the observed rates are 2.8, 
5.1, 2.2, 0.85 and 7.2 for S Per, PZ Cas, VY CMa, NML 
Cyg, and VX Sgr, respectively. By considering the median 
of these numbers, the van Loon rates are ~ 2.8 times too 
large on average. For comparison, the ratios between de 
Jager rates and observed rates for the same stars are 0.20, 
0.78, 0.45, 0.15 and 2.9. The median is 0.45 so that, the de 
Jager rates are ~ 2.2 times too small on average. In conclu- 
sion, the de Jager and van Loon rates are almost equally 
accurate (or inaccurate), but one (de Jager) underestimates 
and the other (van Loon) overestimates the observed rates. 
The de Jager formalism is closer to the observed rates for 
non-dusty RSGs, although it is poor in the case of fj, Cep 
(off by an order of magnitude) . 

4.4. The behaviour of the de Jager rate with effective 
temperature 

Figure 6 shows in particular the temperature dependance 
for the van Loon and de Jager prescriptions. It is surprising 
that for log L ^ 4.5, the M of de Jager gets smaller for de- 
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Table 6. Red supergiants with an observed mass- loss rate 
derived from observations of the gas (M gas ). Column 5 is 
the de Jager rate, and column 6 is the rate with the van 
Loon recipe, which is applicable only for dusty stars (rates 
in parentheses should not be considered). Mass-loss rates 
are given in units of 10 -6 M© yr _1 . 



creasing T e gf. We cannot explain this behaviour in view of 
the data concerning the 15 RSGs that dJ88 used when de- 
riving their formula (these data are listed in Appendix A). 
In particular, VY CMa has the largest mass-loss rate of 
their small sample (2.4 x 10~ 4 M yr _1 ) and its temper- 
ature adopted by dJ88 is low, T c g — 2840 K. In addition, 
the weight assigned to this star by dJ88 is considerably 
larger (60) than the weight of any other RSG in the sample 
(weights ~ 1-4). Yet, the de Jager rate decreases for cooler 
temperatures. A possible explanation of this peculiar T e g 
behaviour is the fact that the fit of the data for the dJ88 
rate was made globally on all stars over the H-R diagram 
and not just "locally" on the RSGs. 

4.5. Is the scatter of M at a given luminosity due to 
effective temperature ? 

An interesting question is whether the deviation of the mea- 
sured M from the dJ88 rate is correlated with effective tem- 
perature. Thanks to the work of LM05, T e g is known for 
the 27 first RSGs of Table 1 (from V589 Cas to TZ Cas). 
T ff is also known for the 8 stars of Table 6. We adopt 
the infrared-based M from Table 1, but prefer when possi- 
ble Mgas from Table 6. Figure 7 (panel a) plots the ratio 
(in log) of measured mass-loss rate to the de Jager rate. 
It shows a slight tendency of M decreasing with increas- 
ing T e ff. The standard deviation is 0.410 when all 31 stars 
are considered and it can be seen that 6 stars (of 31) have 
ratios larger than 4 or smaller than 1/4. Changing the de 
Jager rate by adding a T~^* factor would slightly reduce 
this scatter. Values of a between 6 and 10 are suitable. We 
show in panel b the resulting situation for (T e ff/3600)~ 10 . 
The standard deviation is 0.367 for all objects included. It 
is 0.296 if /i Cep and V396 Cen are omitted, and all other 
objects are within a factor of 4 of agreement. V396 Cen 
has a poor quality 60 ^m flux, coded F in the IRAS point 
source catalog, and presents an anormally low /60//25 flux 
ratio. Therefore its position in Fig. 7 is plausibly explained 
by an underestimated 60 /im flux. 
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Fig. 7. Panel a: Discrepancies between measured mass- 
loss rates (MLR) and mass-loss rates predicted by the de 
Jager prescription, plotted versus T e ff, for Galactic red su- 
pergiants. The horizontal dotted lines indicate discrepan- 
cies by factors 4 and 1/4. Panel b: same as panel (a) when 
the de Jager rate is multiplied by (T e ff/3600)~ 
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4.6. The RSGs in the Magellanic Clouds and the 
dependence on metallicity 

We have seen that the mass-loss rates in the LMC from 
Groenewegen et al. do not follow fully the de Jager rate. 
For L between 63000 and 160 000, there is reasonable agree- 
ment, but at larger L , the de Jager rate is off (too large) by 
a factor of ~ 8. Recently, Bonanos et al. (2010) have also 
considered the mass-loss rates of RSGs in the LMC and 
SMC. Their sample of RSGs (~ 100 sources in the LMC, 
34 in the SMC) is taken from a list of luminous objects 
having optical spectral classification and photometry from 
0.3 to 24 /im obtained during the Spitzer SAGE survey of 
the Magellanic Clouds (Meixner et al. 2006). Bonanos et 
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Fig. 8. Mass-loss rates versus luminosity for Magellanic 
Clouds. The (median) data of Bonanos et al. (2010) are 
represented by histograms (see text). The solid and dotted 
lines show the dJ88 rates for 3500 and 4000 K scaled for 
various metallicity, as indicated in the panels. 



al. (2010) estimate mass-loss with a M versus the K — [24] 
relation derived from model dusty envelopes. 

From the data plotted in their Fig. 17, which show im- 
portant scatter at a given luminosity, we derive median 
mass- loss rates in several bins in L. The 5 bins for the LMC 
have respectively 15, 24, 28, 25 and 5 data points. The 4 
bins for the SMC have 6, 10, 11 and 7 data points. The 
resulting M versus L relations are shown in Fig. 8 (his- 
tograms). Both the SMC and LMC data suggest that M 
increases with luminosity if logi is less than ~ 5.3, as al- 
ready concluded by Bonanos et al. (2010). It can be noted 
that the LMC bin with the highest luminosity suggests that 
M could saturate, but the low number of points makes it 
difficult to be sure of that (see these 5 points at the highest 
luminosity in Fig. 17 of Bonanos et al. 2010). However, we 
note that this is quite consistent with our Fig. 5 and its 
high-L, low-M group of stars, as mentioned above. 

In order to investigate how the RSG mass-loss rates 
depend on metallicity, we can compare the median rates 
obtained above with the dJ88 prescription scaled by various 
Z values. Here we adopt Z(LMC) = Z Q /2, and Z(SMC) = 
Zq/5. Figure 8 shows the de Jager lines for warm effective 
temperatures which are suitable for the Clouds. It shows 
that, for the SMC, the observed rates are between the de 
Jager rate scaled by Z 5 and Z, and a Z - 7 scaling would 
fit the data quite well. As for the LMC, the situation is 
less clear because the metallicity decrease is only a factor 2 
with respect to solar. If we omit the most luminous bin, a 
Z - 7 scaling is not unreasonable, but a Z 5 is also possible. 



5. Conclusions 

The de Jager prescription used in evolution codes for RSGs 
was established in 1988 when the mass-loss rates of stars 
located over all the H-R diagram were fitted as a function 
of L and T c ff ■ The fit was not specifically made for this type 
of stars, and the goal of this work was to compare the de 
Jager rate to estimates of M based on data acquired after 
1988. Our main conclusions are the following: 

1. When the Jura's formula is adopted to get an obser- 
vational estimate of M based on the IRAS infrared excess 
at 60 /im and supposing a gas-to-dust ratio of 200 (corre- 
sponding to complete condensation of refractory material) , 
it is found that the Galactic RSGs are satisfactorily dis- 
tributed around the de Jager line. The deviations from de 
Jager rate is in general less than a factor of 4. 

2. Only eight Galactic RSGs have their circumstellar gas 
detected, i.e. through detection of Hi or CO emission, or 
other methods permitting to estimate M from observation 
of the gas. For these stars, M agrees reasonably well with 
the de Jager rate. The most discrepant cases are found at 
high luminosity, where the de Jager prescription underesti- 
mates M of NML Cyg by a factor of 7, but overestimates 
M for n Cep by a factor of 10 (see Table 3). 

3. The Spitzer data from Bonanos et al. (2010) show that 
the mass-loss rates of RSGs in the Small Magellanic Cloud 
agree on average with the de Jager rate scaled for metallic- 
ity by {Z/Z Q f- 7 . 

4. A similar behaviour is less clear with RSGs of the Large 
Magellanic Cloud. For luminosities smaller than 1.6 x 10 5 
L Q , the median mass- loss rates roughly follow the de Jager 
prescription and its increasing trend with luminosity. These 
measured rates are also roughly compatible with the Z - 7 
scaling found for the SMC. However, for higher luminosi- 
ties, one finds that the mass-loss rates level off at a rate 
of ~ 3 x 10~ 6 M yr _1 and are smaller by a factor of 
<~ 8 than the (unsealed) dJ88 rate (with the exception of 
WOH-G64). 

5. Finally, although the dJ88 prescription presents some 
odd aspects that can be criticized (e.g., its behaviour with 
T e g), it passes rather satisfactorily several tests (listed 1,2,3 
and part of 4 above). Therefore, we think it more appro- 
priate to keep it unchanged in the stellar evolution models 
but to apply a Z 7 factor for metallicity dependence. The 
LMC mass-loss saturation issue needs to be clarified. We 
hope to further contribute to that subject in the future. 
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Appendix A: List of red supergiants considered by 
de Jager et al. (1988) 

The list of RSGs considered by dJ88 are in Table Al. The 
data from dJ88 are: the stellar effective temperature Ti 
(K), the luminosity L\ (L Q ) and the mass- loss rate M\ 
in M Q yr _1 . The columns T 2 , L 2 and M 2 give the same 
quantities as derived in this work. Many mass-loss rates 
used by dJ88 (Mi) came from the exploitation and analy- 
sis of infrared flux measurements between 3.5 and 11.5 /tm 
(Gehrz & Woolf 1971). The mass-loss rates of Hagen (1982) 
are derived from flux measurements at 20 and 25 ^m. The 
mass-loss rates of Sanner (1976) are based on high resolu- 
tion spectral profiles of strong resonance lines. For objects 
with Mi based on several individual determinations (from 
4 to 9 determinations for /j, Cep and a Ori, respectively), 
see the details in de Jager et al. (1988). 

Appendix B: Photometric data of the studied 
RSGs 

Table Bl presents the magnitudes in the UBVIJHK filters 
and the IRAS fluxes at 12, 25, and 60 /im that were used 
to calculate the flux averaged wavelength A m , and the 
luminosities L. In case that a magnitude in the above 
filters was not available, we used the following relations 
which were established from the sample data: 

[/« 0.62 + 0.78 {B-V) 

I » V - 0.80 (V - J) 

K - 0.27 - 0.119 (K - [12]) 
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Table A.l. List and data for the RSGs considered by dJ88. 



Ubject 


'/ ■ 
Ji 


J-2 


logLi 


log 1/2 


!og(MiJ 


10g(iW2) 


Note 


TV Gem 


3560 


3700 


4.75 


4.92 


-5.92 


-6.46 


GW 


BU Gem 


3440 


3800 


5.15 


4.93 


-5.68 


-6.23 


GW 


a Ori 


3500 


3810 


4.87 


4.75 


-6.02 


-5.82 


(1) 


a Sco 


3540 


3550 


4.74 


4.85 


-6.08 


-6.00 


(1) 


RWCyg 


3365 


3600 


5.20 


5.16 


-5.30 


-5.50 


GW 


VV Cep 


3175 




5.13 




-6.89 




GW 


CE Tau 


3400 




4.55 




-6.62 




(4) 


u Cep 


3400 


3700 


4.55 


5.61 


-5.78 


-5.30 


(1) 


Rf Over 


O LOO 


OO 1 o 


o.o^l 




^ 1 ^ 




Vjr vv yz, J 


BI Cyg 


3155 


3575 


5.38 


5.09 


-5.15 


-5.33 


GW 


VY CMa 


2840 


3430 


5.76 


5.47 


-3.62 


-4.40 


(1) 


EV Car 


2930 




5.74 


5.83 


-6.00 


-4.90 


(3) 


S Per 


2810 




5.66 


4.93 


-4.57 


-5.16 


GW 


UY Set 


3000 




5.44 




-5.22 




GW 


AH Sco 


2800 




5.62 




-6.00 




(3) 


NML Cyg 


2485 


3309 


5.80 


5.50 


-4.18 


-3.85 


(1) 



Notes: GW indicates that Mi was taken by dJ88 from Gehrz and Woolf (1971). (1) indicates that Mi is the average of several 
determinations obtained with different methods. (2) BC Cyg has no IRAS fluxes at 25 and 60 /im and was not included in this 
study. (3) Mi from Hagen (1982). (4) Mi from Sanner (1976) 



/bo « 0.20 f 25 
/ioo ~ 0.33 f 60 



Appendix C: Wind expansion velocities 

Figure C.l shows the wind expansion velocities of several 
RSGs plotted versus luminosity. The Galactic RSG data 
are from Table 2 of this paper. The LMC RSG data are 
from Marshall et al. (2004) (their Table 1 for luminosities 
and their Table 2 for expansion velocities). The error bars 
represent typical uncertainties on luminosity, ± 30% for 
the Galaxy and ±20 % for the LMC. The straight lines 
have the slope predicted by Habing et al. (1994) for dust- 
driven winds and were used in the Jura's formula for red 
supergiants lacking a wind velocity measurement. 



LMC RSGs 





1 4.5 5 5.5 

Log (Luminosity) 

Fig. C.l. Wind expansion velocity, in kms -1 , as a func- 
tion of luminosity (in Lq), for Galactic RSGs (lower 
panel), and LMC RSGs (upper panel). The straight line 
for Galactic RSGs corresponds to V = 20 (L/10 5 ) - 35 in 
kms -1 . The straight line for LMC RSGs corresponds to 
V = 14 (L/10 5 ) ' 35 . 



13 



N. Mauron and E. Josselin: The mass-loss rates of red supergiants 

Table B.l. Adopted photometry for calculating L and A m . The UBVIJHKL arc in magnitudes, and the fluxes / are in 
Janskys 



Obiert 


u 


B 


v 


I 


j 




H 




K 


L 


J 12 


for 


f 

J 60 


BD+60335 = V589 Cas 


13.94 


11.49 


9.15 


5.03 


4.06 


3 


.01 


2 


.49 


2.16 


12.74 


8.67 


3.61 


BD+56512 = BU Per 


14.35 


11.70 


9.23 


4.80 


3.68 


2 


.68 


2 


.19 


1.85 


44.97 


31.29 


5.23 


HD 14469 = SU Per 


12.04 


9.80 


7.63 


3.82 


2.82 


1 


.93 


1 


.46 


1.13 


48.71 


30.66 


6.87 


HD 14488 = RS Per 


12.90 


10.62 


8.35 


4.33 


3.05 


2 


.11 


1 


.56 


1.22 


74.44 


47.82 


9.93 


HD 14528 = S Per 


14.51 


11.88 


9.23 


4.48 


2.95 


1 


.84 


1 


.12 


0.55 


339.40 


233.20 


40.59 


HD 14826 = V441 Per 


13.02 


10.56 


8.24 


4.30 


3.47 


2 


.47 


2 


.04 


1.60 


18.53 


13.25 


3.54 


HD 936979 = YZ Per 


13.14 


10.55 


8.20 


4.34 


3.26 


2 


.30 


1 


.91 


1.64 


38.85 


26.12 


5.28 


W Per 


14.80 


12.24 


9.62 


4.75 


3.09 


2 


.00 


1 


.57 


1.51 


90.58 


78.86 


14.87 


BD+57647 


14.93 


12.26 


9.52 


4.76 


3.83 


2 


.71 


2 


.11 


1.56 


39.04 


26.45 


6.47 


HD 37536 = NO Aur 


10.52 


8.30 


6.21 


3.01 


2.12 


1 


.13 





.97 


0.51 


43.49 


22.89 


5.12 


a Ori 


4.36 


2.29 


0.40 


-2.47 


-2.93 


-3 


.73 


-4 


.01 


-4.43 


4680.00 


1740.00 


299.00 


HD 42475 = TV Gem 


10.58 


8.81 


6.56 


3.17 


2.20 


1 


.20 





.94 


0.54 


96.08 


41.16 


6.06 


HD 42543 = BU Gem 


11.11 


8.63 


6.39 


3.07 


2.23 


1 


.28 





.98 


0.62 


78.20 


47.59 


10.50 


CD-314916 = V384 Pup 


13.37 


11.07 


8.91 


5.43 


4.61 


3 


.49 


3 


.14 


2.83 


10.85 


7.34 


2.76 


HD 90382 = CK Car 


12.00 


9.66 


7.45 


3.55 


2.63 


1 


.71 


1 


.36 


1.03 


113.40 


70.10 


13.98 


HD 97671 = V602 Car 


13.50 


10.91 


8.39 


3.61 


2.49 


1 


.53 





.99 


0.09 


174.90 


85.13 


12.40 


V396 Cen 


12.30 


10.00 


7.85 


3.58 


2.58 


1 


.56 


1 


.20 


0.79 


53.30 


39.81 


4.98 


KW Sgr 


14.81 


12.50 


9.66 


4.23 


3.12 


1 


.99 


1 


.54 


0.91 


250.10 


148.30 


18.39 


iliJ ooUyUo4 — IN xx VUl 


1 K A 1 
10.41 


10/11 
lz.41 


y.ou 


A AG 
4.4y 


O.ZO 


O 

z 


1 A 
. 14 


1 




i nn 
l.UU 


i nA Qn 
lUu.oU 


Do. oU 


1 OQ 


BD+364025 = BI Cyg 


14.83 


12.34 


9.33 


3.64 


2.35 


1 


.15 





.58 


0.02 


334.60 


244.90 


51.23 


KY Cyg 


17.44 


14.63 


11.14 


4.01 


2.35 





.85 





.17 


-0.46 


510.60 


329.60 


50.74 


BD+394208 = RW Cyg 


13.43 


11.01 


8.13 


3.38 


2.11 





.94 





.45 


-0.08 


298.40 


189.80 


60.69 


H Cep 


8.88 


6.43 


4.17 


0.31 


-0.52 


-1 


.30 


-1 


.65 


-2.09 


1296.00 


607.70 


127.00 


Case 75 = V 354 Cep 


16.95 


13.85 


10.67 


5.20 


3.92 


2 


.57 


1 


.87 


1.39 


77.67 


46.56 


8.01 


Case 78 = V355 Cep 


15.47 


13.06 


10.76 


5.78 


4.61 


3 


.29 


2 


.68 


2.28 


19.04 


14.47 


3.27 


BD +602613 = PZ Cas 


12.84 


11.48 


8.90 


3.90 


2.42 


1 


.53 


1 


.02 


0.39 


373.00 


398.20 


96.48 


BD +602634 = TZ Cas 


14.17 


11.68 


9.17 


4.62 


3.34 


2 


.31 


1 


.88 


1.37 


75.07 


51.70 


9.47 


EV Car 


12.43 


10.09 


7.89 


3.21 


2.11 


1 


.21 





.79 


0.32 


265.40 


164.00 


25.87 


HS Cas 


14.52 


12.33 


9.82 


5.13 


3.98 


2 


.88 


2 


.46 


2.05 


23.47 


16.41 


3.51 


XX Per 


11.72 


10.36 


8.26 


4.25 


3.06 


2 


.21 


1 


.83 


1.19 


70.76 


28.42 


4.23 


KK rer 


12.36 


9.9/ 


1.(6 


a on 


3.00 


o 
z 


1 A 

.14 


1 


.68 


1. (0 


19.3/ 


10.8/ 


o no 
A. LA 


AD Per 


12.66 


10.13 


7.85 


4.30 


3.38 


2 


.48 


1 


.94 


1.65 


21.28 


13.90 


2.85 


PR Per 


12.80 


10.14 


7.84 


4.60 


3.56 


2 


.68 


2 


.25 


2.14 


12.79 


9.39 


2.37 


CP Cas 


14.67 


12.08 


9.37 


4.76 


3.53 


2 


.38 


1 


.95 


1.51 


25.73 


18.47 


4.45 


VY CMa 


12.01 


10.19 


7.95 


3.10 


1.96 





.43 


-0 


.73 


-2.68 


9919.00 


6651.00 


1453.00 


a Sco 


4.08 


2.78 


0.96 


-1.92 


-2.55 


-3 


.47 


-3 


.70 


-4.04 


3198.00 


689.90 


115.50 


VX Sgr 


11.72 


9.41 


6.52 


2.11 


1.23 





.13 


-0 


.50 


-1.61 


2738.00 


1385.00 


262.70 


Case 49 


17.33 


14.26 


11.12 


5.41 


4.07 


2 


.74 


2 


.15 


2.01 


17.86 


11.01 


6.58 


U Lac 


12.50 


11.04 


8.70 


4.30 


2.90 


2 


.13 


1 


.57 


0.77 


124.00 


61.49 


9.04 
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